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Abstract. We present tiie results of a detailed temporal analysis of the bright BL Lac object Mrk 421 using the 
three available long timing mode observations by the EPIC PN camera. This detector mode is characterized by 
its long life time and is largely free of photon pile-up problems. The source was found in different intensity and 
variability states differing by up to more than a factor of three in count rate. A time resolved cross correlation 
analysis between the soft and hard energy bands revealed that the characteristics of the correlated emission, with 
lags of both signs, change on time scales of a few 10^ seconds. Individual spectra, resolved on time scales of ~ 100 s, 
can be quite well fitted by a broken power law and we find significant spectral variations on time scales as short as 
~ 500 — 1000 sec. Both the hard and the soft band spectral indices show a non-linear correlation with the source 
flux. A simple power law model of the form F ex flux~" with ahard ~ 0.13 and asoft ~ 0.22 describes rather well 
the observed trend of decreasing F values with increasing flux, which appear to "saturate" at the same limiting 
value of F ~ 1.8 at the highest flux levels. A comparison of the observed light curves with numerical results from 
relativistic hydrodynamic computer simulations of the currently favored shock-in-jet models indicates that any 
determination of the jet's physical parameters from 'simple' emission models must be regarded with caution: at 
any time we are seeing the emission from several emission regions distinct in space and time, which are connected 
by the complex hydrodynamic evolution of the non-uniform jet. 

Key words. BL Lacertae objects — individual: Mrk 421; X-rays: galaxies — Radiation mechanisms: non-thermal. 
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1. Introduction (z = 0.031) object was observed by almost all previous X- 

^ ray missions and shows remarkable X-ray variability corre- 

BL Lacs are radio-loud AGN dominated by relativis- lated with strong activity at TeV energies (e.g., Takahashi 

tic jets seen at small angles to the line of sight et al. 1996, Maraschi et al. 1999). Structure function and 

d ' ( |Urry fc Padovani 1995| ). They are characterized by large power density spectrum analyses indicate a roll-over with 

luminosities, large and variable polarization, and strong a time scale of the order of 1 day or longer (Cagnoni et 

variability. A substantial part of the information we pos- al. 2001, Kataoka et al. 2001) which seems to be the time 

sess about these objects is obtained from the analysis scale of the successive flare events. On shorter time scales 

of temporal variations of the emission and the combined less power in the variability is found with a steep slope 

spectral and temporal information can greatly constrain of the power density spectrum ^ /^(^"'^^ (Tanihata 2002, 

the jet physics. Time scales are related to the crossing Brinkmann et al. 2003). 
time of the emission regions which depend on wavelength 
and/or the time scales of the relevant physical processes. 

The measured lags between the light curves at different The smooth featureless continuum spectra of Mrk 421 

energies as well as spectral changes during intensity varia- from the radio to the X-ray band imply that the emis- 

tions allow us to probe the physics of particle acceleration sion in these energy bands is due to relativistic particles 

and radiation in the jet. radiating via the Synchrotron process; the hard X-rays 

Mrk 421 is the brightest BL Lac object at X-ray and and 7-rays are likely to be produced by inverse Compton 

UV wavelengths and it is the first extragalactic source dis- scattering of synchrotron photons by the same electrons 

covered at TeV energies (Punch et al. 1992). The nearby (Ulrich et al. 1997 and references therein). Broken power 

law models provided better fits to the ASCA 0.6—7 keV 

Send offprint requests to: W. Brinkmann; e-mail: data than a simple power law (Takahashi et al. 1996) with 

wpb@mpe.mpg.de a break energy of 1.5 keV, and the change of the power 
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law index at the break point is AF ~ 0.5, i.e., the spec- 
trum is steeper at higher energies. 

With the wider energy range of BcppoSAX it be- 
came clear that the simple models are not adequate de- 
scriptions of the downward curved Synchrotron spectra 
(Fossati ct al. 2000b) and continuously curved shapes had 
to be employed (Inoue & Takahara 1996, Tavecchio et al. 
1998). The Synchrotron peak energy varied between 0.4— 
1 keV, the spectral index at an energy of 5 kcV between 
2.5 < F < 3.2. Both quantities are correlated with the X- 
ray flux: with increasing flux the synchrotron peak shifts 
to higher energies and the spectrum gets flatter at higher 
energies. 

The emission of the soft X-rays is generally well corre- 
lated with that of the hard X-rays and lags it by 3—4 kscc 
(Takahashi et al. 1996, 2000, Zhang et al. 1999, Malizia 
et al. 2000, Kataoka et al. 2000, Fossati et al. 2000a, Cui 
2004), however, significant lags of both signs were detected 
from several flares (Tanihata et al. 2001, Takahashi et al. 
2000). 

These results were obtained from data with relatively 
low signal-to-noise ratio, integrated over wide time inter- 
vals (typically one satellite orbit) and by analyzing data 
from prominent flares with time scales of a day. Recently, 
XMM — Newton with its high sensitivity, spectral resolv- 
ing power and broad energy band provided uninterrupted 
data with high temporal and spectral resolution. In an 
analysis of the early XMM — Newton data (Brinkmann 
et al. 2001) for the flrst time the evolution of intensity 
variations could be resolved on time scales of ~ 100 s. 
Temporal variations by a factor of three at highest X- 
ray energies were accompanied by complex spectral vari- 
ations with only a small time lag of t = 265j^}Q2 s be- 
tween the hard and soft photons. Scmbay ct al. (2002) 
confirmed these short lags in an analysis of a larger set of 
XAdM — Newton observations of Mrk 421 and they show 
that the source exhibits a rather complex and irregular 
variability pattern - both temporarily and spectrally. From 
XMM - Newton observations of PKS 2155-304 Edclson 
et al. (2001) suggest that previous claims of soft lags with 
time scales of ^ hours might be an artifact of the peri- 
odic interruptions of the low-Earth orbits of the satellites 
every ^1.6 hours. This claim was questioned by Zhang 
et al. (2004) who show that, although periodic gaps intro- 
duce larger uncertainties than present in evenly sampled 
data, lags on time scales of hours cannot be the result of 
periodic gaps. 

In an extended analysis of all available 
XMM - Newton PN observations of Mrk 421 
Brinkmann et al. (2003) split up the individual light 
curves into shorter sub-intervals in which the flux 
appeared to show some similar specific behavior. The 
cross-correlation analysis of these individual soft and hard 
band light curves demonstrated that they are in general 
well correlated at zero lags, but sometimes the hard band 
variations lead the soft band variations by typically ~ 5 
min, in two cases the soft band lead the hard band varia- 
tions. The delays appear to be correlated with a harder 



spectrum during higher intensities. Recently, Ravasio 
et al. (2004) used three XMM — Newton observations 
in Nov/Dec. 2002 when Mrk 421 was highly variable. 
During two large flares the source showed strong spectral 
evolution and temporal lags between the diflFerent energy 
band variations of both signs, with lags of up to more 
than 1000 sees. 

In this paper we will present a detailed analysis of the 
three currently available XMM data sets on Mrk 421 in 
which the PN camera was operated in timing mode. We 
re-processed the data and used the most recent response 
functions. In previous observations of Mrk 421 even the 
high time resolution of the PN in small window mode 
sometimes proved to be insufficient to prevent photon pile 
up when the source was in a high intensity state. This 
led to a reduction of the signal to noise ratio for the spec- 
tral data and to frequent unpredictable large observational 
gaps in the light curves, due to the extremely high work 
load of the on-board electronics (FIFO overflows) and the 
limits of the telemetry rate (counting mode). Timing mode 
data have a very high photon collecting efficiency of 99.5% 
of the frame time and do not suffer from pile-up problems 
which would, with other detector modes, lead to a sub- 
stantial loss of data for the extremely high count rates of 
Mrk 421 . 

Two of the data sets, that of orbit 84 as well as the 
short flare in orbit 546 have been presented previously 
(Brinkmann et al. 2003 and Ravasio et al. 2004, respec- 
tively), however, only in a more limited analysis and with- 
out special emphasis on their short time scale behavior. 
The third dataset, from orbit 807, when the source is 
extremely bright, has not been published before. In our 
analysis we take full advantage of the excellent statistical 
quality of the timing mode data sets and, in particular, 
we will for the first time analyse the data on time inter- 
vals considerably shorter than the individual observations 
as the source exhibits changes in its spectral properties on 
time scales down to a few thousand seconds. We will study 
the spectral slope variability on time scales as short as ^ 
1000 sec and how that correlates with the flux variations; 
we search the cross correlation function over very short 
periods and discuss how it changes with time as the flux 
varies. 

We will first present the data and demonstrate the 
spectral complexity of Mrk 421 in broad band spectral 
fits. We will discuss the 'standard' cross-correlation anal- 
ysis of the data which is commonly used to investigate the 
cross-links between the hard and soft band light curves 
and will then introduce a sliding-window technique to 
study the cross-correlations on shorter time scales. In Sect. 
4 we will present a detailed analysis of the source proper- 
ties using temporally resolved spectral analyses. We then 
discuss these results in the framework of the shock-in-jet 
model with some emphasis on numerical simulations of the 
hydrodynamic evolution of non-uniform relativistic jets. 
Final conclusions will be given in Sect. 6. 
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Orbit 


Observing date 




PN mode 


Filter 


Duration 


Live Time 


count rate 


source count rate 




(UT) 








[ks] 


[ks] 


on chip 


[0.6-10 keV]* 


084 


May 25, 2000: 03:53 - 


10:11 


Timing 


Thick 


22.4 


16.78 


373.3 


221.6±0.1 


546 


Dec 01, 2002: 23:18 - 


18:44'!' 


Timing 


Medium 


71.2 


59.00 


303.9 


176.9±1.3 


807 


May 06, 2004: 02:38 - 


21:11 


Timing 


Medium 


66.7 


25.63 


841.3 


491.7±0.2 



: Next day 

^: Background subtracted source count rate for the spectral fit 



2. The XMM-Newton observations 

Mrk 421 was observed repeatedly during the XMM- 
Newton mission, mostly as a calibration source for the 
RGS. During three of these observations the PN camera 
was operated in timing mode, once with a thick filter (or- 
bit 84), the other times with medium filter. This mode has 
the advantage that even at highest source count rates, pile 
up, and thus spectral distortions, can be neglected. The 
details of the observations are given in Table ^ where we 
first list the orbit in which the observation took place, the 
date of the observation, the instrument settings and the 
duration of the observation. Due to the high count rate 
from Mrk 421 the detector falls frequently into counting 
mode, which introduces gaps of typically 21 s duration into 
the data stream. This happens when the read-out buffer 
is 75% full, and the rate at which this occurs depends 
critically on the work load of the detector. For example, 
in orbit 807, with the extremely high count rate on the 
chip, about every 20 sec of data are followed by this ~ 21 
sec gap; the observation consists of ^ 1800 of these 'ele- 
mentary observation intervals' (EOIs). In orbit 546, with 
its much lower count rate, each EOI consists of ^ 100 
sec data, followed by the gap. Correspondingly, the Live 
Time of the detector (column 6) is rather short for orbit 
807 and highest for orbit 546. Column 7 gives the average 
count rate recorded on the chip when no spatial nor qual- 
ity selection for the photons were made, while in column 8 
we list the screened, background subtracted source count 
rates for the total spectral fit in the 0.6—10 keV band. 

All data have been reprocessed using XMMSAS ver- 
sion 6.1. For the spectral fits the latest response matrices 
were used, for which adjustments of the effective areas 
around the Au edges were made as well as corrections for 
the thick filter near the Al edge. We selected photons with 
Pattern < 4 (i.e. singles and doubles) and quality flag — 
0. For the source data we selected the photons from rows 
29 < RAWX < 45 for our analysis, as background we used 
the photons from rows 2 < RAWX < 18 (for a description 
of the instrument and the detector modes see Ehle et al. 
2001). The energy range for all analyses was taken to be 
600 eV < S < 10 keV as below 600 eV the single-double 
ratio is not constant and thus gives unreliable results in 
the spectral analyses. The background count rates in the 
whole energy band are generally less than ~ 1.5% of the 
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Fig. 1. Background subtracted, 0.6—10 keV PN light 
curves of Mrk 421 with time binning of 100 s. Times are 
counted from the beginning of the actual exposure. The 
curves are labeled by the orbit of the observation. 



source count rates, and even at energies ^ 10 ke'V they 
amount to ^ 10% of the source rates. A few short back- 
ground flares and the data at the end ( ^ 60 ksec) of the 
observation of orbit 807 were discarded from the spectral 
analysis presented in Sect. 4. 



2.1. Orbit integrated properties 

To demonstrate the magnitudes of the source variability 
we present in Fig.nthe 0.6—10 keV 100 s bin light curves 
for the three observations with the above photon selection 
criteria. The counting errors are smaller than the sym- 
bol sizes. For comparison with previously published light 
curves please note that the count rates are dominated by 
photons at low energies. For example, extending the en- 
ergy range down to 500 eV would increase the count rate 
by about 15%, going down to 200 eV nearly doubles the 
number of counts. The orbit 84 observation was performed 
with the thick filter; correcting these data to the medium 
thick filter would result in a factor of typically ^ 1.2 
higher count rates. 
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Fig. 2. Broken power law model fit to the 0.6—10 keV data 
during the three observations. The two top figures show the fit 
and the data to model ratio for orbit 84, the middle panel the 
ratio for orbit 546, the bottom panel that for orbit 807. 

The light curves are characterized by long term trends 
lasting for time scales comparable to or longer than the 
length of the observations and, superimposed on this, of 
frequent flare like events of relatively regular shape with 
characteristic time scales of the order of a few thousand 
seconds. Interestingly, larger variations occur when the 
source is in a higher state: during orbit 546 the maximum 
to minimum intensity variations are ^ ± 15%, in orbit 84 
of - ±20% and in orbit 807 even > ± 40%. 

2.2. Spectral analysis 

As mentioned above the spectral behavior of Mrk 421 is 
rather complex and with increasing data quality more so- 
phisticated spectral models had to be used for the fits. To 
show the integral spectral behavior of the source during 
the three observations we present in Fig.|21the fits to the 
total data of each orbit to a broken power law. The plot- 
ted residuals are heavily rebinned. The deviations between 
model and data are only of a few percent. The numerical 
results are given in Table 2. The absorption has been fixed 
to the galactic value; keeping Nh as a free parameter re- 
sults in slightly lower values of x^d ^^'^ slightly higher 
column densities. Interestingly, the slope changes at the 
break energy are rather small which might indicate that 
we see the source near the Synchrotron peak of its SED. 



Table 2. Results for the broken power law fits assuming 
galactic absorption (Nh = 1-66 x 10^° cm~^) 



Orbit Tsoft Ebrk Thard Norm x?cd/dof 
(keV) (1) 

084 2.257±0.001 4.00±0.17 2.354±0.009 0.182 1.22/1584 
546 2.400±0.009 3.37±0.42 2.620±0.006 0.108 1.69/1671 
807 1.969±0.001 4.63±0.19 2.027±0.006 0.316 1.58/1900 



(1): Normalization at 1 keV in ph/keV/cm^/s. 



Fitting a continuously curved model (Fossati et al. 2000b) 
to the data results in fits of very similar Xredi however the 
residuals indicate that at higher energies the curvature of 
the spectrum is over-estimated. 

The Xrod given in Table 2 show that the fits are hardly 
acceptable which is to a great extent due to the extremely 
high photon statistics in the soft energy band (for exam- 
ple, there are 13.2 million photons in the orbit 807 spec- 
trum). There is still some uncertainty in the calibration, 
visible around ^ 2 keV (the Au edge and, possibly, some 
Si feature at slightly lower energies) and a general con- 
vex curved low energy part. Whether this slight extra low 
energy curvature is in the range of the calibration uncer- 
tainties or a 'source intrinsic' spectral property cannot be 
answered yet. The residuals around 2 keV further indi- 
cate that these structures are not strictly stationary dur- 
ing the three observations taken over four years. They are 
likely not source specific but they could be small secular 
changes of the detector characteristics over these years. 
(A gain-fit over a limited energy range around 2 keV re- 
duces the residuals strongly) . Secondly, the spectral slopes 
vary with the intensity variations causing some disper- 
sion of the data with respect to the average model. The 
time resolved fits of a broken power law model to shorter 
data segments (see Sect. 4) mostly result in Xrcd ^ l-Oj 
thus a broken power law provides an adequate description 
of these Mrk 421 spectra. The general impression from 
the above total fits is that the spectrum gets flatter with 
increasing flux and that the break energy shifts towards 
higher values. 

3. Cross-Correlation Analysis 

In order to investigate the cross-links between the hard 
and soft band light curves we estimated their Cross- 
Correlation Function (CCF), using light curves with a 
10-s bin size. If a data segment is followed by the typi- 
cal 21 s long FIFO gap, this is filled with a running mean. 
We originally chose to work with this small bin size be- 
cause the signal is rather strong and, secondly, previ- 
ous cross-correlation analyses between soft and hard band 
light curves have indicated that, if inter-band delays exist, 
they are very small (Brinkmann et al. 2003). The follow- 
ing temporal and spectral analyses show, however, that 
noticeable spectral changes occur on longer time scales of 
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Fig. 3. CCF plots for the orbit 84 (top), orbit 546 (middle) 
and orbit 807 (bottom) observation. The 10 s binned soft and 
hard band light curves for the total observations were used; 
the error bars are smaller than the symbol sizes. The curves 
are labeled by the orbit of the observation. 

a few hundred seconds. Since the results from the previous 
section show that the energy spectrum exhibits a signif- 
icant spectral slope change at energies between 3 and 4 
keV (see Table 2), we chose the count rates in the 0.6—2 
keV and 4—10 keV bands as representative of the soft and 
hard band spectral components of the source. This choice 
of energies further ensures that there is no 'spill-over' of 
flux between the two bands. 

We calculated the sample Cross-Correlation Function, 
CCF{k), as follows: 

nm?(h\ - T.ti^soft{t) - Xsoft){Xhard{t + k) - Xhard) 

k = 0,±At,...,±{N -l)At. 

The summation goes from t = At to {N — k)At for 
fe > (At ~ 10 s, and is the total number of points 
in the light curve) and the sum is divided by the number 
of pairs included, i.e. N{k). For negative lags fc < the 
summation has to be done over x^oftit+ \ k\) Sind :x.hardit)- 
The variances in the above equation are the source vari- 
ances, i.e., after correction for the experimental variance. 
Significant correlation at positive lags means that the soft 
band variations are leading the hard band. 

First, we computed the CCF (up to lags ±3 hrs) us- 
ing the whole length (i.e. total) soft and hard band light 
curves. The resulting CCFs are shown in Fig.|21(the CCF 
for orbit 84 is similar to Fig. 13 of Brinkmann et al. (2003), 



while Ravasio et al. (2004) present in Fig. 9 the CCF of 
only the 'flare' at t 42 ksec in orbit 546). Although all 
three CCFs peak at zero lag, there are clearly definite dif- 
ferences in shape between them. The ranges for the values 
of the cross-correlation coefficients are different. Further, 
in orbit 84 the auto-correlation function of the soft band is 
very similar to the CCF; in the other two orbits the CCF 
is skewed with respect to the ACF in the sense that the 
soft/hard bands are less correlated at negative lags, but 
more correlated at positive lags 1 hour. 

This behavior, a peak at lag zero and a skewed CCF, 
has been observed many times. Ravasio et al. (2004) 
showed that such a shape can be be reproduced by as- 
suming flares that have the same form, a linear rise 
and exponential decay, but with different amplitudes and 
rise/decay time scales in the two energy bands. This expla- 
nation certainly works in the case of well sampled isolated 
flares. Here, we do not focus on single events but rather 
use the whole light curve to estimate the CCFs. The ob- 
served fight curves are characterized by many "events" 
rising or decaying on different time scales, suggesting that 
more than one emission region is in operation in the sys- 
tem. The emission mechanism in one of them may oper- 
ate on time scales equal to/longer than the duration of 
the light curves. It causes the larger amplitude variations 
which happen in phase in both energy bands, hence the 
zero lags when we integrate over the whole length of the 
observation. On the other hand, emission from other re- 
gions may operate on shorter time scales, causing changes 
of the emission pattern during the observation on a "char- 
acteristic time" scale of just a few ksec. The physical 
meaning of this characteristic time scale is not obvious. 
Obscrvationally, in the present data it could correspond 
to the duration of the flare- like events of relatively regular 
shape in the light curves which is equal to a few thou- 
sand ksec. It is possible that this time scale corresponds 
to the size of the individual emission regions in the source. 
Weaker, shorter (i.e. on time scales of seconds) variability 
might remain undetected in the measurement noise. On 
the other hand, there are also longer time scales in the 
system as, for example, the large scale variations during 
orbit 807 demonstrate. We now investigate whether the 
correlation properties of the emission in the two energy 
bands do evolve on periods comparable to this 'few thou- 
sand second' time scale. 

Long integration intervals were required in the past 
to obtain an acceptable signal to noise ratio with the 
limited photon statistics of previous instruments. The 
largely increased sensitivity of XMM — Newton allows 
the analysis of much shorter data streams. Brinkmann et 
al. (2003) split the individual Mrk 421 observations into 
two or three sub- intervals with similar variability patterns 
and flux levels and found differing CCFs for these sub- 
intervals. 

Here we will employ a 'sliding window' technique and 
calculate the CCFs for shorter data intervals, which start 
at different times of the observed light curve and range 
over a restricted time interval. Thus, the above defined 
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Fig. 4. Sliding window CCFs for orbit 546. The amplitude 
of the CCF is color coded, with lags plotted in the vertical 
direction. The top figure was calculated with an individual data 
length of 2000 s, the middle with 6000 s, and the bottom figure 
with 20000 s long data streams. The bottom panel shows the 
color bar used for the numerical values of the CCFs for the 
above figure and for Fig. ^ 



CCF(k) will then be replaced by a CCF(k; T, £), i.e., the 
cross-correlation coefficient at a lag 'k' is calculated for a 
data stream with length C which starts at the time T in 
the light curve. 

In Fig. ^ we show a two dimensional representation 
of these sliding window CCFs, based on light curves with 
10 s binnings. The vertical scale represents the lags (sim- 
ilar to the X-axis of Fig. 3), the color coding indicates 
the amplitude of the cross-correlation coefficient, with the 
color code given in the lowest panel of the figure. Time is 
along the x-axis. In the top figure the length of the indi- 



0.2B r 




5.0x103 l.QxlO* 1.5x104 2.0x10^ 
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Fig. 5. Standard deviations (Jp^ccf) of the sliding window 
CCFs for orbit 546 for different data lengths £. The diamonds 
are the measured values, the three large symbols indicate the 
window lengths used for the CCFs in Fig. 0] The solid line 
gives the mean and the dashed lines the standard deviation of 
the re-sampled CCFs (see text). 

vidual data sets is £ = 2000 s (thus the longest part of 
the light curve is covered), for the middle it is 6000 s and 
for the bottom figure 20000 s. Clearly, in the first figure, 
CCF(k; T,2000), the cross-correlation coefficient is very 
noisy, the CCF is dominated by local fluctuations and the 
data length C seems to be shorter than the typical variabil- 
ity time scale. The middle panel {C = 6000 s) shows vary- 
ing lags of both signs with different amplitudes as well as 
periods where there is hardly any inter-band correlation. 
This length seems to be in the range of the variability time 
scale of the source and it further demonstrates that the 
physical conditions of the emission regions change on that 
time scale. In the bottom figure the integration length 
of £ = 20000 s over the data stream is obviously much 
longer than the variability time scale. Local variations of 
the sign and amplitude of the lags are smoothed out and 
produce cross-correlations with lags around zero, with an 
asymmetric amplitude distribution as found in Fig. 13 

The above discussion shows that there is an optimal 
window length C which reveals best the time evolution 
of the CCF structure. In order to determine this length 
we regard the values ('response') of the 'local' cross cor- 
relation functions, CCF(k), in all sliding windows for lags 
-2000 s < k < +2000 s as a probability distribution func- 
tion (PDF), P{ccf), of different response values ranging 
from -1 to 1. By determining the statistical properties of 
these PDFs as a function of C we can identify the opti- 
mized size of C. If the sliding window is very small, only 
small correlations between the bands can be identified, be- 
cause the CCF is dominated by the random fluctuations 
of the experimental noise. This results in a relatively nar- 
row PDF peaking at P{ccf) « 0. If £ is very large, one 
is no longer sensitive to the temporal fluctuation of the 
CCF, which may vary with appearing or disappearing flare 
events. Only an overall 'mean cross correlation' between 
the bands is measured. The resulting PDF is also narrow 
but with a maximum at high values accounting for the 
obvious fact that there is an overall correlation between 
the energy bands. However, if the sliding window has an 
intermediate 'optimal' length one is sensitive to the varia- 
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tions of the CCF, if present. The resulting P[ccf) becomes 
broad with CCF-values covering a large part of the whole 
possible range and no pronounced maximum. To quantify 
this anticipated behavior of the PDFs we calculated their 
standard deviation crp(j,c/) and plot them as a function 
of the sliding window sizes, L (Fig. |3J). The shape of the 
curve shows a broad maximum with the highest value at a 
window length C ^ 7.2 ksec. At window sizes greater than 
~ 15 ksec the standard deviations gradually decrease. 

To verify that the observed behavior is of physical na- 
ture and not caused by statistical effects, namely less sam- 
pling for larger window sizes, we performed the following 
validation, which is inspired by ideas from bootstrapping 
(Efron 1993): For each window size we randomly selected 
as many CCFs as we have for the largest window size, 
and use them as the basis for the calculation of the stan- 
dard deviation. This procedure enables us to determine 
a kind of bootstrap-error of the calculated api^^^f^ by re- 
peating the re-sampling procedure many times and cal- 
culating statistical means and standard deviations of the 
measured quantity sigma. The results are shown in Fig. 
|S1 as the three lines (mean and standard deviations). It 
is obvious that we obtain the same results as previously, 
when the number of windows changed with £. 

Therefore, the highest variations in the CCF(k; T, £) 
at a window length C ^ 7.2 ksec truly indicate that this 
value matches the typical time scale in the light curve and 
that cross correlation functions calculated over the total 
light curve (i.e. ^ 20 ksec) provide only little information. 

For the two other orbits a similar behavior was found 
with slightly different values for the maximal values (orbit 
84: peak at 5.2 ksec; orbit 807: peak at 10 ksec). In Fig. 
|S|we show the CCFs for all three orbits calculated with 
these three optimized window sizes, for orbits 084 (top), 
orbit 546 (middle) and orbit 807. Again the amplitude of 
the CCF is color coded, with lags plotted in the vertical 
direction. Below the individual CCFs we show the cor- 
responding light curves in the total (0.6—10 keV) energy 
band. The x-axis is the time from the beginning of the 
observation which, for the CCFs, corresponds to the start 
time of the CCF calculation. The light curves end at the 
time given by the difference of the observing time and the 
length of the sliding window. 

The plot clearly demonstrates the complex behavior 
of the energy resolved light curves and the considerable 
uncertainties interpreting them. Let us note some pecu- 
liarities which are hard to accommodate in 'standard' sce- 
narios of the emission from the BL Lac jets: 
Orbit 084: Generally, there are always nearly zero lags 
between the two energy bands. The correlation coeffi- 
cient tends to zero from times t ~ 2—3 ksec, at the small 
bump in the light curve, and at around 8— 9. 5 ksec, at 
the continuously increasing part of the light curve. As 
the data length C of orbit 084 is taken to be 5.2 ksec, 
the CCF at ^ 2 ksec, where it changes from highly cor- 
related to weakly correlated, is calculated from data be- 
tween 2— 7. 2 ksec; that where it changes back between ^ 
3— 8. 2 ksec. Therefore, the data segments that can cause 
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Fig. 6. Sliding window CCFs for orbits 084 (top), orbit 546 
(middle) and orbit 807 (bottom) calculated with the optimized 
sliding window lengths (see text). The amplitude of the CCF is 
color coded, with lags plotted in the vertical direction. Plotted 
below the individual CCFs are the corresponding light curves 
in the total (0.6—10 keV) energy band. The x-axis is the time 
from the beginning of the observation which, for the CCFs, 
corresponds to the start time of the data window for the cal- 
culations. 
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those changes are either between 2— 3ksec or between 
7.2— 8.2ksec, i.e., at the beginning or at the end of the 
data sets. The hght curve does not provide a clear an- 
swer to this question; however, the time resolved spec- 
troscopy (see next section. Fig. IT)) clearly points to a short 
hard flare centered at t ^ 3.5ksec. The final parts of the 
CCF at times ^ 10^ s, which cover the intensity drop by 
about 25% at t ^ ISksec and the subsequent recovery 
of the count rate indicate that these flux variations have 
occurred at nearly zero lags. 

Orbit 546: A small positive lag at the beginning, followed 
by zero lags during the following plateau phase. The in- 
tensity drop at t = 20ksec results in a strong negative 
lag, while the increase and the subsequent decrease at t 
XL 38—45 ksec results in positive lags. The final parts 
of this decrease and the small substructure at ~ 51 ksec 
seem to have occurred nearly uncorrelated at the two en- 
ergy bands. Note that the absolute variability amplitude 
during this observation is relatively small. 
Orbit 807 : The orbit where the source is brightest and ex- 
hibits the largest variability. Correlated at near-zero lags 
at the beginning and then, around the turning point of 
the flux increase, an extended period of very low correla- 
tion. The large variability region a.t t ^ 30 ksec seems to 
have occurred in both energy bands without strong time 
lags. However, the broadness of the correlation function 
at these later periods could indicate that the time scales 
of the smaller intensity fluctuations are not resolved. 

These investigations lead us to the conclusion that, 
quite generally, an analysis of light curves in diflerent 
energy bands might be inconclusive; in particular, when 
rapid and strong intensity variations cannot be resolved 
temporarily, for example, due to a low signal to noise ra- 
tio. The above CCF results demonstrate the complexity of 
the observed variations of the light curves in different en- 
ergy bands but do not provide further information on the 
spectral behavior. An intensity decrease/increase is once 
accompanied by soft/hard lags (respectively) and once oc- 
curs nearly achromatically. A more detailed understanding 
of the variability behavior of the source can be achieved 
with a study of the variations of the spectral slopes on 
time scales as short as possible, which we present below. 

4. Time-resolved spectroscopy 

The observed count rates in the PN camera are high 
enough to study the energy spectrum of Mrk 421 over 
short periods of time in order to investigate in detail its 
spectral variability properties. To this end, we extracted 
energy spectra over the short observation intervals be- 
tween the gaps that are introduced when the detector falls 
into counting mode (see Sect. 2 above). The typical dura- 
tion of these intervals are ^ 100 s, 80s and 20 s for the ob- 
servations of orbit 84, 546 and 807, respectively. To avoid 
insufhcient signal for extremely short EOIs we had to add 
in some cases subsequent data intervals to get the above- 
mentioned 'typical' data stretches. Thus, for orbit 84 we 
fltted 149 individual data segments, 476 for orbit 546 and 



456 for orbit 807. The number of photons for the indi- 
vidual flts in the 0.6 — 10 keV band were between 15000 
and 30000, which is large enough to perform a meaningful 
model fitting to each spectrum individually. We used the 
same response matrices as for the total energy spectrum 
(see Sect. 2.2) and again fitted a broken power law model. 

As a first test of the physical relevance of our results we 
investigated the best fit normalization values. We found 
an excellent correlation between these values and the ob- 
served count rate in the 0.6—10 keV energy band indicat- 
ing that the observed flux variations are mainly caused by 
variations in the power law normalization (as opposed to 
spectral slope variations). 

In most cases, the best fitting break energy value 
turned out to be in the range between 1.5 — 5 keV, i.e. 
similar to the values we found for the fits to the total en- 
ergy spectra of the source. In some cases though, the full 
band energy spectrum could be quite well fitted by a single 
power law, and thus either the hard or soft spectral slopes 
i^hard Or Tsoft, respectively) were determined with sig- 
nificant uncertainty. Furthermore, in most cases the best 
fitting Thard and Tgoft values are not different by more 
than ~ 1 — 2a. This is due to the fact that the uncertainty 
associated with the hard band spectral slope is often quite 
large which is a direct result of the fact that most of the 
photons in the individual spectra are in the soft energy 
band. However, as our aim is to investigate the intrinsic 
temporal variations of the spectral "shape" of the source, 
the somewhat larger uncertainties of the individual best 
fitting slope values are not so important. What is impor- 
tant is to investigate whether the best fitting values are 
systematically different between successive data stretches, 
and if so, in what way. 

Due to their large uncertainties and the existence of a 
few outliers, the light curves of the best fitting hard and 
soft slope values are noisy. For that reason we binned the 
respective light curves by calculating the weighted mean 
of 4 consecutive slope estimates in the case of orbits 84 
and 807 observations and of 12 values for orbit 546. As a 
result, we are able to study the spectral slope variations 
of the source on time scales as short as ^ 500 s and ^ 
700 s in the case of orbits 84 and 807, and ^ 180 s for 
orbit 546. The bottom panels in Figs. [TH^l show plots of 
the binned Fhard.so/t light curves. The soft and hard band 
count rate light curves are plotted in the top panels of the 
same figures. 

These figures show clearly that apart from the strong, 
fast flux variations, Mrk 421 also exhibits significant spec- 
tral variations on time scales as short as ^ 500 — 1000 sec. 
To the best of our knowledge, this is the first time that 
these quantitative spectral variations on time scales of a 
few hundred seconds have been studied along with the flux 
variations. We observe significant slope variations with a 
max-to-min variability amplitude of the order of ~ 7—10% 
and ~ 4 — 6% in the hard and soft slope light curves, re- 
spectively, for orbits 84 and 546. Apart from the largest 
amplitude flux variations, orbit 807, when the source is at 
its brightest state, shows the largest amplitude spectral 
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Time (sec) 

Fig. 7. Background subtracted, 0.6-2 and 4-10 keV PN 
light curves of Mrk 421 in orbit 84 (open circles and 
filled squares, respectively), in bins of size 200 sec (top 
panel). The light curves are normalized to their mean. The 
corresponding Fhord and Tsoft curves (filled squares and 
open circles, respectively) are plotted in the bottom panel 
(binned as described in the text). The shaded boxes iden- 
tify those parts of the observation where the T vs source 
flux relation is studied in more detail (see Sect. 4.1). 



variations as well: the max-to-min amplitude is ~ 30% 
and ~ 20% for the Thard and Tgoft hght curves, respec- 
tively. In all cases, the hard band slope variations are of 
larger amplitude when compared with the soft band slope 
variations. Furthermore, both of them are of smaller am- 
plitude than the amplitude of the flux variations in the 
respective bands. 

Apart from the significant, fast spectral variations, the 
bottom panels of Figs. [7|- El show that, most of the time, 
the hard and soft band slopes are not the same. This result 
justifies a posteriori the use of the broken power law model 
fit to the individual, time resolved spectra of the source. 
Although the uncertainties of the slope values from the fits 
to each individual spectrum are not negligible the avail- 
ability of many spectra, and the use of a broken power 
law model, reveals clearly that the soft and hard band 
slopes differ in a systematic way. In most cases, the slope 
of the photon spectrum at hard energies is steeper than 
the soft band slope. There are cases though when the pho- 
ton spectrum is well described by a single slope power law 
(for example the period between 20—30 ksec after the start 
of the observation in orbit 807) and even cases when the 
hard band slope is flatter than the soft band slope (for 
example the period after the first 30 ksec in orbit 807). 

The light curves of the fitted break energies show sig- 
nificant variations which are neither correlated with the 
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Fig. 8. Same as Fig. [7\ but for the observation during 
orbit 546. The hard band and soft band count rate light 
curves are in bins of size 1000 s and 500 s, respectively. 
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Fig. 9. Same as Fig. [7\ but for the observation during 
orbit 807. The count rate light curves are in bins of size 40 
s. 

best fitting slopes nor with the sources' flux states. This 
may well be representative of the intrinsic behavior of the 
source. On the other hand, as discussed in Sect. 2.2, the 
X-ray spectrum of Mrk 421 may show a "continuous" 
curvature which the broken power law model can only ap- 
proximate to some level. Thus, the break energy may not 
correspond to an intrinsic characteristic of the source, but 
may simply "divide" the spectrum into the two regions 
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that show the largest difference in "curvature". For this 
reason, we do not study the break energy variations here- 
after. Instead, we concentrate on the study of the Thard 
and ^soft variations, having in mind that these slopes sim- 
ply give a measure of the spectral curvature towards the 
soft and hard end of the X-ray spectrum. 



4.1. The relation between spectral and flux variations 

"Hardness (or softness) ratio vs count rate" plots have 
been used in the past mainly to study the correlation be- 
tween the flux and spectral variations in this source and 
other BL Lacs (see for example Brinkmann et al. (2003) 
and Ravasio et al. (2004) for a recent use of these plots 
with XMM — Newton data). In some cases, the evolu- 
tion of the spectral slope itself as a function of the source 
flux has also been investigated. For example, Takahashi 
et al. (1996) studied the spectral slope variations of the 
source based on a day-long ASCA observation, Fossati et 
al. (2000b) used the best continuously curved model fits 
to study the spectral evolution on time scales of ^ 10 
ksec, based on BeppoSAX observations, while Ravasio et 
al. (2004) have used simple power law fits to study the 
spectral variations during a few particular periods, chosen 
from two XMM — Newton light curves. 

We believe that the analysis presented in the previous 
section provides a further improvement in these studies. 
Due to the high sensitivity of XMM — Newton we can 
probe spectral variations on time scales as short as ^ 1 
ksec, based on model fitting results rather than on hard- 
ness ratios. By fitting a broken power law we are able to 
investigate separately the relation between the slopes of 
the soft and hard bands of the spectrum with the count 
rates in the respective energy bands. Also, we will study 
the cross-correlation between the soft and hard band slope 
light curves. 

Figs. EHni show that the hard and the soft band slope 
variations are correlated to a large degree. In general, both 
parts of the spectrum steepen or flatten at the same time, 
although not always at the same rate. There are cases 
when the hard and soft band variations are somehow "dis- 
connected". For example, in the period around 7 ksec, be- 
tween 11 — 13 ksec and at 18 ksec during orbit 84, we 
observe fast, large amplitude Thard variations, which are 
absent in the T^oft time series. Similar events are observed 
in orbit 546, and at 7—11 ksec in orbit 807. Not surprising, 
during the same periods, the peak in the CCF between the 
soft and hard band light curves is decreased (Sect. 3). 

In order to investigate this anti-correlation further, we 
plot in Fig. 1101 the slopes Fhard and ^soft as a function 
of the soft and hard band count rates, using all the data 
shown in Figs. [7HH1 The top and middle panel in this 
figure show clearly the strong anti-correlation between the 
spectral and flux variations: both the hard and soft band 
spectra become harder as the flux decreases. This is a 
well known behavior for this source and other BL Lacs 
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Fig. 10. The Thard vs 4 — 10 keV count rate and the 
Fso/t vs 0.6 — 2 keV relations (top and middle panels, re- 
spectively). The bottom panel shows the ratio T hard soft 
plotted as a function of the 0.6 — 10 keV count rate. Data 
points from Figs.CUHl and^lare shown with open circles, 
open squares and filled triangles, respectively. The solid 
curves describe broad trends as discussed in the text. For 
the orbit 84 observation, which was performed with a thick 
filter, the 0.6 — 2, 4 — 10 and 0.6 — 10 keV count rates are 
multiplied by factors of 1.22, 1.01, and 1.19, respectively. 
These are the count rate ratios expected for the thick and 
medium thick filters in the case of a power law spectrum 
and Nh similar to that of Mrk 421 . 



(c.f. Plan 2002, and references therein; Brinkmann et al. 
2003). 

The fiux - spectral slope relation is obviously not lin- 
ear. Both the hard and the soft band slopes appear to 
"saturate" to the same limiting value of F ~ 1.8 at the 
highest flux levels. This has already been observed in the 
past, see e.g. Ravasio et al. (2004). However, the large 
number of data points in Fig. 1101 suggest a more compli- 
cated picture. For example, although most points in the 
spectral slope vs flux plots do follow a power law like rela- 
tion (shown with the solid lines in the two upper plots) , in 
other cases a linear relation appears to describe the data 
better. It seems possible that there are various well de- 
fined paths (i.e. "spectral states") in the "F vs fiux" plane 
that the source follows at different times. 

The bottom panel in Fig. ^| shows a plot of the ra- 
tio Thard/Y soft vcrsus thc total (i.e. 0.6 — 10 keV) source 
count rate. We find that the slope ratio does not remain 
constant, but rather decreases with increasing source fiux. 
A power law model (shown as a solid line in this panel) 
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Fig. 11. The average cross-correlation function between 
Thard and 4—10 keV count rate (top panel), the r^o/t 
and 0.6 — 2 keV count rate (middle panel), and between 
the two spectral slopes Thard and Tsoft (bottom panel), 
using data from all observations. The CCFs are defined 
such that any positive delays would imply a delay of the 
spectral with respect to the flux variations (in the upper 
two panels), and a delay of the Thard with respect to the 
Tsoft variations, in the bottom panel. 

appears to be broadly consistent with the decrease of the 
slope ratio as the source flux increases, although not all 
points follow the same trend. 
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Fig. 12. Orbit 84: the Thard vs 4 - 10 keV count rate and 
the Tsoft vs 0.6 — 2 keV relations (top and middle panels, 
respectively) for the "PI" and "P2" parts (identifled by 
the shaded boxes in Fig. Cj). The bottom panel shows the 
ratio Thard/Tsoft as a function of the 0.6 — 10 keV count 
rates. Filled circles show the data points during the flux 
rise and open circles during the decaying phase of each 
event. The dashed line indicates the linear path that the 
rising (or the decaying) phase points follow. For clarity, in 
cases where the two phase points follow different paths, or 
deflne loop-like structures, the respective points are con- 
nected with solid lines. 



4.2. Temporal correlations of the spectral variations 

To investigate whether the spectral and flux variations 
happen simultaneously, we computed the cross-correlation 
function (CCF) between the spectral slope and count rate 
light curves shown in Figs. [7\-^ We used the "Discrete 
Correlation Method" of Edelson & Krohk (1988), with a 
lag size of Afc = 800 sec, 1000 sec, and 500 sec in the 
cases of orbits 84, 546 and 807, respectively. We also com- 
puted the DCF between the T^ard and Tsoft hght curves, 
to examine the simultaneity of the different band slope 
variations. Note that in this case we did not estimate the 
CCF at lag zero, because the slopes were obtained from 
model fitting to the same spectrum and any correlated, 
systematic biases in the estimation of Tsoft and T^ard can 
artificially increase the CCF value at this lag. 

The cross correlation functions of all three observa- 
tions look very similar. For that reason, in order to in- 
crease the signal-to-noise, we combined the three individ- 
ual Thard/Tsoft VS flux and the Thard vs Tsoft CCFs into 
one file, respectively, and estimated average CCFs in each 
case. The average Thard vs fiux, Tsoft vs flux and Thard 



vs Tsoft CCFs are shown in Fig. ^2 The highest peak is 
at lag ~ in all cases, suggesting that any delays between 
the slope and flux variations, or between the soft and hard 
band slope variations, should be shorter than ~ 1000 sec. 
The CCFs look fairly symmetric around zero lag suggest 
that both the hard and soft band spectra and flux vari- 
ations happen almost simultaneously, to the limit of the 
current temporal resolution. 

4.3. The shortest time scales 

Although the results from Sect. 4.1 show that, in gen- 
eral, the spectrum hardens as the source fiux increases, 
the r vs fiux behavior of the source is quite complicated 
during individual "events" . For that reason, we now focus 
our study of the spectral variability behavior of the source 
during the periods which are marked with shaded boxes in 
Figs.lTHHl These correspond to well defined, single "flare- 
like" events, where both the flux rising and decaying parts 
have been adequately sampled (with the exception of the 
P2 event in orbit 546, where it is not clear whether the 



12 



W. Brinkmann et al.: XMM— Newton timing analysis of Mrk 421 




no 
1.1 

1.08 
I 1.06 
' 1.04 

1.02 



130 120 130 

0.6-2 keV Count Rate (c/s) 




150 140 
0.6-10 keV Count Rate (c/s) 



Fig. 13. Same as Fig. [El but for the parts "PI" and "P2" 
in orbit 546. 
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Fig. 14. Same as Fig. [HI but for the parts "PI" and "P2" 
in orbit 807. 



flux has started decreasing or not). Figs. ll2H14lshow the 
^hard vs 4 — 10 keV count rate (top), the Tgoft vs 0.6 — 4 
keV (middle) and the Thard/Tsoft vs 0.6 - 10 keV (bot- 
tom panels) during these events. Filled and open circles 



correspond to data points during the flux rise and decay 
phases, respectively. 

The results from the analysis presented in this section 
are complementary to those presented in the past (e.g. 
Takahashi et al. 1996, Fossati et al. 2000b, Ravasio et al. 
2004). However, they also provide some new insight into 
the source's behavior. By a) studying all flaring events 
that we could identify in the XMAI — Newton light 
curves, by b) using spectral slope estimates in two en- 
ergy bands on time scales as short as ~ 1 ksec, and by c) 
combining spectral with timing information (from Sect. 3) 
we demonstrate that the source does not operate in the 
same way at all times. 

Our results in the case of the PI flare during orbit 546 
are in agreement with the results of Ravasio et al. (2004). 
Although our model fitting results show clearly that the 
the soft and hard band parts of the spectrum have differ- 
ent slopes, both Thard and r^o/t follow a clockwise loop 
pattern during the evolution of the flare. In fact, even 
the Thard/^ soft ratio follows a clockwise loop as the flare 
evolves. As the flux rises, both the soft and hard band 
spectra flatten, but at different rates: ^soft flattens faster. 
When the flux decreases the spectra become steeper, but 
P/iard steepens faster than Tgoft, and the slope ratio de- 
creases for a while. Finally, it increases again as Thard does 
not vary appreciably any more while T soft still increases. 
The CCF results suggest hard lags during the same period 
(see Sect. 3), a result that is consistent with the spectral 
evolution behavior of the source: variations are propagated 
from the softer to the harder energy bands during the ini- 
tial flux rising phase, and then the hard band flux decays 
faster than the soft band flux. 

However, the source does not show the same behavior 
at all times. For example, during the PI flare in orbit 807 
the Thard variations follow the same path during the flux 
rise and decay phase (shown with the dashed line in the 
top left panel of Fig. I14|l . No loop patterns are observed 
in the Thard/T soft ratio either; the slope ratio remains 
roughly constant, implying that the hard and soft band 
spectra evolve at the same rate. The CCF results imply no 
delays between the variations observed in the two bands. 
A similar behavior is also observed during the P2 flare in 
orbit 84. Both Thard and T soft evolve in the same way 
during the rising and decaying phases (dashed lines in the 
top and middle right panels in Fig. 1121) . Consequently, the 
slope ratio remains roughly constant as the flare evolves. 
In agreement with the observed spectral variability behav- 
ior, the CCF analysis shows no delays between the hard 
and soft band during this event. 

The PI flare in orbit 84 is a case where obvious vari- 
ations in Thard are not evident in the T soft time series. 
The dashed line in the top left panel in Fig. 1121 shows the 
Thard variations as the flux increases (filled circles). Some 
of the Thard values during the fiux decay phase do fol- 
low the same path, but most of the Thard values during 
the prolonged " fiare-plateau" phase between 11 — 13 ksec 
(Fig.EJ) are clearly not consistent with it. Furthermore, the 
Tsoft values during fiux decay lie systematically above the 



W. Brinkmann et al.: XMM— Newton timing analysis of Mrk 421 



13 



Tso/t values in the rising phase suggesting the presence of 
an anti-clockwise loop pattern. As a result of these two 
effects, the slope ratio also follows a loop pattern, which 
evolves in an anti-clockwise direction. The P2 flare in or- 
bit 807 is another case of a flare where the Tfiard and T^oft 
variations do not show the same behavior: the former fol- 
lows a clockwise loop pattern (top right panel in Fig. [T^i , 
while the ^soft evolution follows the same linear path in 
both the rising and decaying flare phases. 

That obvious case where the hard and soft band slopes 
are almost entirely "disconnected" is the P2 event in or- 
bit 546. Initially, T^ard steepens (AF ^ —0.05) and Tg^ft 
flattens (AF ~ 0.05) although the flux does not change 
appreciably. Then, both slopes remain constant during a 
~ 15% and ^ 8%, rapid hard and soft flux increase, re- 
spectively. Towards the end of the observation, Tgoft re- 
mains constant, while Thard flattens slightly (right pan- 
els in Fig. I13|l . However, the slope ratio changes imply 
that these apparently different variations are in fact "con- 
nected" in some way. A well defined loop pattern appears 
in the slope ratio evolution plot which implies that the 
variations are propagated from the soft to the hard band. 
This result is also supported by the strong, hard lags that 
are observed in the same period (see Sect. 3). 

5. Discussion 

In this work we presented a spectral and timing analysis 
of the three calibration phase XMM — Newton observa- 
tions of Mrk 421 performed in timing mode. The analysis 
strongly relied on that mode as, for example, the enor- 
mously high count rates of ~ 800 cts/s in orbit 807 would 
have lead to a substantial photon pile up - even in the 
extremely short frame time of ~ 5.7msec of the FN small 
window mode and thus to a large reduction of photons 
available for the analysis. 

Our average results confirm and, in a way, expand the 
findings of previous XMM — Newton investigations of 
Mrk 421 (Brinkmann et al. 2001, Sembay et al. 2002, 
Brinkmann et al. 2003, and Ravasio et al. 2004) and ear- 
lier extended studies with ASCA and BeppoSAX. In the 
following we summarize our main findings which might 
help to understand the complex behavior of the source. 

5.1. The CCF results 

In general, the observed soft and hard band variations are 
well correlated in all three observations. However, using a 
sliding window cross-correlation analysis, we could show 
that the cross links between the soft and the hard bands 
are really complicated and not easily understood. There 
exists a 'characteristic' time scale (of the order of a few 
ksec) on which the CCF appears to change "continuously" 
and, depending on the length of the observing window and 
the actual activity state of the source, we find periods with 
positive or negative or no lags, but also periods of weak 
correlations between the soft and hard energy bands: 



— Feriods of weak correlations appear when variations 
in the hard band flux or hard band spectra slope are 
absent in the soft band. 

— Strong correlation with no obvious lags is observed 
throughout orbit 807, although the fiux does go up and 
down. In fact, this orbit shows the highest amplitude 
flux variations! 

— In cases where we observe the soft leading the hard 
band variations (e.g. PI and P2 in orbit 546) or the 
hard leading the soft (e.g. at ~ 23 ksec in orbit 546) 
the delays are just a few minutes in agreement with 
the lags already noticed by Brinkmann et al. (2003). 

— The analysis of orbit 84 shows a much more compli- 
cated picture than anticipated from Brinkmann ct al. 
(2003) where the observation was split in only two 
parts. Here we find the hard fiux leading the soft at 
the beginning, then loss of correlation, then the soft 
fiux leading, then the hard leading briefly, and then 
a loss of correlation in the prolonged "flare plateau" 
phase, where the hard band variations are absent in 
the soft band. 

There is, certainly, a continuous range of variability 
time scales, as seen for example in orbit 807, and up to a 
day or longer (Kataoka et al. 2001). However, the above- 
mentioned 'characteristic' time scale of correlated emis- 
sion, i.e., the time scale over which the correlation between 
hard and soft band emission is found to exist, to change 
sign, or to disappear, might be representative for the size 
of the individual emission regions. With a bulk Lorentz 
factor of the jet of F ^ 10 this would correspond to an 
emission region of a fewxlO^^cm. 

5.2. The relation between spectral and flux variations 

The cross-correlation analysis demonstrates, on the other 
hand, the limitations of this method as only little informa- 
tion about the actual spectral variations can be extracted. 
We therefore used the high signal to noise ratio to perform 
proper spectral fits to the data down to time scales of ^ 
20 sees. Mrk 421 shows significant spectral variations on 
time scales as short as ^ 500 — 1000 sec and most of the 
time, the hard and soft band slopes are not the same. 
Generally, the slope of the photon spectrum at hard ener- 
gies is steeper than the soft band slope but there are cases 
when the hard band slope is flatter than the soft band 
slope 

The spectral behavior follows the well-known 'the 
brighter the harder' trend. Both, the hard and the soft 
band spectral indices show a non-linear correlation with 
the source flux, and there is an indication that the source 
may populate different spectral states at different times 
in the "F vs ffux" plane. In general, as the flux rises, both 
the soft and hard band spectra flatten, but 7 soft flattens 
faster. Then the spectra become steeper as the fiux de- 
creases, but now Vhard stcepcns faster than T soft- On long 
time scales ^ soft and Thard vary in phase and they vary 
in phase with fiux. 
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On shorter time scales, during the evolution of some 
flare like events, we observe "hysteresis" phenomena. 
^softj^hard and their ratios follow loop- like paths in the 
spectral vs flux plots in the clockwise (PI in orbit 546) 
and anti-clockwise (PI in orbit 84) direction. In other 
cases though, both evolve on a straight line during rise 
and decay of the flux (P2 in orbit 84), we observe a flare 
where only one follows a loop (P2 in 807), and there is an- 
other case where only the slope ratios follow a well deflned 
structure (P2 in orbit 546). 



5.3. Implications to physical models 

The numerous instances where the source does not fol- 
low the "average" standard pattern show that any deter- 
mination of the physical parameters of the emission re- 
gion from simple 'homogeneous', 'one-zone', etc., models 
might be questionable. There are cases, for example PI 
and P2 in orb 546, which could be "flares" that corre- 
spond to emission from a single, perturbed region of the 
jet. They show hard lags, so the whole events are gov- 
erned by the "acceleration" time scale of the relativistic 
electrons. Even the soft lags that we see during the flux 
drop before PI, could indicate evolution governed by the 
cooling time scale. Interestingly, these "clear" events show 
up when the source is at its lowest flux state. It could be 
that in these cases we have mainly only one region emit- 
ting at any time. In general, the complex spectral behavior 
indicates that the emission geometry is not that simple. 
At a time tobs we are not only seeing the emission from a 
certain position R(t) in the relativistically moving jet, but 
as well the earlier emission from all positions ~ R - c (5t 
down the jet. In the framework of the commonly accepted 
shock-in-jet model this emission is determined by the his- 
tory of the colliding shells (blobs) - and, perhaps, by the 
earlier (or later) emission of other, independent shells. 

Even in the framework of the basic, 
two— colliding— shells model (Spada et al. 2001), non- 
linearities are expected. Most of the emission is, certainly, 
generated when the second, faster shell hits the flrst, 
slower shell. Particles are accelerated, magnetic flelds are 
amplified and part of the available energy of this shocked 
region is radiated. However, the flrst shell travels as well 
with supersonic speeds along the jet, a shock structure 
has developed at its front and a reverse shock travels 
back through the shell and, after some time, runs into 
the shocked emission region of the two colliding shells. 
Depending on the initial (and boundary) conditions of 
the system this might lead to noticeable changes of the 
physical state of the emission region and to 'unexpected' 
changes of the emission characteristics. 

These effects are illustrated in the following flgures, 
obtained from numerical relativistic hydrodynamic simu- 
lations of two colliding shells with Lorcntz factors Fi = 5 
and T2 = 7.5 (Mimica et al. 2004a, 2004b). Fig. [Tsl shows 
the density distribution of the shells at five different times 
in the jet frame as seen by a distant observer. The origin of 
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Fig. 15. Density (in arbitrary units) of the colliding shells 
at different evolution times (in the jet frame) of the system 
as seen by a distant observer. The observation time (with 
arbitrary origin) runs from right to left. The plot shows at 
what times the observer would 'see' the various structures 
approaching him with relativistic velocities. 
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Fig. 16. Normalized X-ray emission from the model of 
Fig.Elin the top panel in the hard (1—10 keV, red curve) 
and the soft (0.1—1 keV) X-ray band. The bottom panel 
shows the evolution of the spectral power law indices in 
the two bands. 



the observer's time scale is arbitrarily chosen. The config- 
uration at T = in the jet frame (black curves) would be 
seen by a distant observer between -10000 < tobs < s as 
two individual shells moving towards him. The situation 
at T ~ 519 ksec in the jet frame (red curve) would flrst 
appear to the observer as the shock structure at the front 
of shell 1 and then, around tobs ^ 2500 sec, as strongly 
rising emission as the two shells have collided. From then 
on, up to about tots ^ 5000 sec, the emission would be 
rising due to the growing density and size of the collision 
region in the jet frame. At tohs ~ 7500 sec, at the evolu- 
tionary time T ~ 852 ksec in the jet frame (yellow curve), 
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yses between the soft and hard energy bands with a sUding 
window technique which revealed that the characteristics 
of the correlated emission changes on time scales of a few 
ksec. With time-resolved spectral fits on time scales of a 
few hundred seconds we could study the spectral evolu- 
tion of the source. We find significant spectral variations 
on time scales as short as ~ 500 — 1000 sec and non-linear 
correlations between between slopes of the fitted broken 
power law model and the fiux of the source. 

The temporal and spectral behavior of the source is 
very complex and during flares various variability pat- 
terns in the soft and hard energy band were observed. 
Correspondingly, it appears hard to deduce uniquely the 
underlying physical parameters for the emission process 
from the observations. We compared the observations with 
relativistic hydrodynamic simulations of the currently fa- 
vored 'shock-in-jet' model for the BL Lac emission (see, 
for example, Spada et al. 2001) which takes into account 
that we are seeing the emission from multiple shocks which 
have either very different physical parameters or that we 
detect the emission from similar shocks at very differ- 
ent states of their evolution, heavily confused by rela- 
tivistic beaming and time dilatation effects. As the ac- 
tual emission characteristics of these numerical models 
depend on a large number of physical parameters and 
poorly known boundary conditions it is a tedious task 
to reproduce numerically the large variety of observed 
light curves. However, BL Lac jets are relatively simple 
physical systems, which might be solely governed by the 
equations of relativistic MHD and special relativity. Thus, 
further observational progress (in particular long continu- 
ous XMM — Newton observations with high spectral and 
temporal resolution) and extended numerical simulations 
might finally lead to a better understanding of the physi- 
cal conditions in these systems. 
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the shock structure from shell 1 has merged with the inter- 
action region of the shells, changing strongly the emission 
properties. Then, at later times, the structure fades away. 

Fig. 1161 shows the normalized light curve of this event 
in the 0.1 — 1 kcV and the 1—10 keV X-ray band in the 
upper panel and the evolution of the spectral indices in 
the lower panel. It must be emphasized that the param- 
eters chosen for this simulation are not at all applicable 
for the physical conditions in Mrk 421. The spectra indi- 
cate that the shocks are too strong and/or the magnetic 
field amplification has been over-estimated; they are just 
presented for illustrative purposes. In particular, the spec- 
tral hardening phase occurs only at the very beginning of 
the rising light curve, then the spectrum starts cooling. 
At the time where the reverse shock from the first shell 
joins the inter-shell emission region, we find spectral and 
intensity changes of the emission which would not fit into 
a 'regular' evolution pattern of the observed flare and the 
cross correlation analysis of this light curve would predict 
lags of changing sign and amplitude. Looking at the above 
presented observations it appears that, for example, parts 
of the hght curve of orbit 807 (Fig. O resemble this be- 
havior: the fiux and spectral index of the soft energy band 
varies smoothly while the flux and power law index in the 
hard band exhibits some 'unexpected' deviations. 

One of the most interesting questions raised by the 
observed emission is the cause of the long-term inten- 
sity changes. In the three observations analyzed here, the 
flux changed by a factor of more than three. Are these 
variations caused by geometrical changes of the viewing 
conditions, like in the lighthouse model by Camenzind & 
Krockenberger (1992)? Or are we seeing periods of en- 
hanced activity, a dramatic increase of the number of in- 
teracting blobs or stronger interactions - for which the 
individual emission patterns are averaged and form an en- 
hanced 'background' emission. There is growing evidence 
that the blazar emission consists of a 'quasi-stationary' 
component which changes only on long time scales and a 
flaring component which gives rise to the short time scale 
intensity and spectral variations. The influence of these 
different components on the determination of the physical 
parameters deduced from the observations has been dis- 
cussed in more detail by Fossati et al. (2000a) and Fossati 
et al. (2000b). In extensive numerical simulations of the 
shock-in-jet scenario, varying the injection rate, the dura- 
tion, speeds and other relevant physical parameters of the 
shells, we could not reproduce the observed variability pat- 
tern of the emission of Mrk 421. Only after adding a sub- 
stantial 'background' emission of ^ 60% to the numer- 
ically estimated flux we did find light curves that looked 
similar to the observed ones. 

6. Conclusions 

We have presented a time resolved analysis of the cur- 
rently available XMM-Newton observations of Mrk 421 
performed with the PN in timing mode. The high signal 
to noise ratio allowed us to perform cross-correlation anal- 
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